Interaction-induced derealization of two particles in a random potential: Scaling 

properties 
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The localization length £2 for coherent propagation of two interacting particles in a ran- 
dom potential is studied using a novel and efficient numerical method. We find that the 
enhancement of £2 over the one-particle localization length £1 satisfies the scaling relation 
£2/^1 — /(^/A^), where u is the interaction strength and the level spacing of a wire 
of length £1. The scaling function / is linear over the investigated parameter range. This 
implies that £2 increases faster with u than previously predicted. We also study a novel 
mapping of the problem to a banded-random-matrix model. 



PACS numbers: 72.15.R, 71.30. 

While much is known about the localization properties 
of one particle moving in a random potential Q, there 
are few secured results about localization in the pres- 
ence of interactions between the particles ||. In view 
of the complexity of the interplay of disorder and inter- 
action, Dorokhov || and, very recently, Shepelyansky 
H approached this problem by studying a simple spe- 
cial case — two interacting particles in a random one- 
dimensional potential — and predicted that the interac- 
tion can lead to a significant derealization of the pair. A 
possible realization of this system are excitons in a dis- 
ordered semiconductor Furthermore, understanding 
the localization properties of two particles in a random 
potential may lead to new insights into the role of inter- 
actions in the Anderson insulator. 

Shepelyansky considered the motion of two particles 
interacting by a short-range interaction in a random po- 
tential. Whenever the two particles are localized far 
apart compared to the one-particle localization length, 
the effect of the interaction is only exponentially small. 
However, an interesting effect occurs when the two par- 
ticles are localized within about one one-particle local- 
ization length of each other. In this case, Shepelyansky 
constructed an approximate mapping of the problem to a 
banded-random-matrix model. Studying this model nu- 
merically, he predicted that independently of the statis- 
tics of the particles and the sign of the interaction there 
is the possibility of coherent propagation of the two par- 
ticles over distances £2 much larger than the one-particle 
localization length £1 M. He found that 
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where u denotes the interaction strength and t is the hop- 
ping matrix element. £1 is measured in units of the lattice 
constant. Subsequently, Imry || has given a Thouless- 
type scaling argument supporting and generalizing this 
result, and Frahm et al. fl have studied the problem 
numerically using a transfer-matrix technique, finding a 



slower increase of £2 with £1 than predicted by Eq. (nl). 
Related results had been found earlier by Dorokhov ||] 
for the propagation of two harmonically bound particles 
in a random potential. 

In this paper, we present a novel and efficient numerical 
technique to compute the two-particle localization length 
£2 directly from a microscopic model. This method al- 
lows us to obtain accurate results over a wide range of 
parameters for both bosons and fermions. Our main re- 
sult is that £2 obeys the scaling relation £2/^1 = f(u/A^), 
where A^ is the single-particle level spacing of a wire of 
length £1. We conjecture, based on our numerical re- 
sults, that the exact scaling function / is linear at the 
center of the band. While our results qualitatively con- 
firm the prediction that a short-range interaction can 
lead to coherent propagation of the pair over distances 
much larger than the one-particle localization length, this 
scaling relation is inconsistent with the original predic- 
tion, Eq. (|l|). It is an important consequence of our 
results that the enhancement of £2 sets in for weaker 
interactions than previously predicted. We also derive 
and study a novel mapping of the problem to a banded- 
random-matrix model. A combination of scaling with 
this banded-random-matrix model suggests the possibil- 
ity that the validity of our principal results extends to 
quasi-one-dimensional wires. 

Our starting point is the Anderson Hamiltonian Hq 
for two spinless particles in a one-dimensional random 
potential with an additional Hubbard-type interaction U, 

H = t {|n, m)(n+l, m\ + \n, m) (n,m+l \ + h.c.} 

n,m 

+ ^\n,m){V n + V m ){n,m\ + U. (2) 



The random site energies V n are drawn uniformly from 
the interval [—W/2,W/2\. The hopping matrix element 
t will be set to unity in the following. We parameter- 
ize the disorder by the one-particle localization length |l| 
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£i = I05(t/W) 2 in the absence of the interaction. For 
bosons, we choose an on-site interaction with matrix el- 
ements (n, m\U\n', to') = u 5 n<n > S m , m >S n . m , for spinless 
fermions a nearest-neighbor interaction with matrix el- 
ements (n,m\U\n',mf) = u[£ n ,ro+i + 5 n ,m-i]Sn,n'S m ,m'- 
In the following our numerical method is described for 
bosons. The extension to fermions is straightforward. 
To study the two-particle localization properties of the 
Hamiltonian (|^) we focus on the matrix elements of the 
two-particle Green function 



g = (e-h - uy 



(3) 



between doubly-occupied sites \n,n). We define the two- 
particle localization length £2 for coherent transport by 
the exponential decrease with distance of these matrix 
elements, 
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This identification is certainly reasonable as long as the 
doubly-occupied sites do not effectively decouple from 
the remaining Hilbert space because of the interaction, 
i.e., as long as the interaction strength u is smaller than 
or of the order of the hopping matrix element t. Our nu- 
merical approach exploits the observation that, for the 
interaction U , a closed equation can be derived for these 
matrix elements. This reduces the dimension of the rele- 
vant Hilbert space from A^ 2 to iV (with iV the number of 
sites) which enables us to study systems with up to 1000 
sites. The Dyson equation for the two-particle Green 
function is 



G — Gq ~\~ GqU G, 



(5) 



where Go = (E—Hq)^ 1 denotes the Green function in the 
absence of the interaction. The interaction can be writ- 
ten as U — uP, where P denotes the projector onto the 
doubly-occupied sites, P\n, to) = 5 ntm \n, to). Hence, mul- 
tiplying the Dyson equation by P on both sides and writ- 
ing U = uP 2 , one obtains a closed equation for the two- 
particle Green function projected onto doubly-occupied 
sites, 



G — Go ~t~ uGqG. 



(6) 



Here we defined G = PGP and G Q = PG P- Solving this 
equation for G one has 



G = 



Go 



1 



l/u - Go 



(7) 



In the site basis, the unperturbed Green function Go is a 
banded matrix whose matrix elements decrease exponen- 
tially with distance on the scale £1 /2. Therefore, we com- 
pute £2 using only the second factor in (|^) from which any 
long-range behavior of G must arise. This is very useful 



for numerical purposes because this factor can be inter- 
preted as the Green function of the "Hamiltonian" Go at 
"energy" l/u. This enables us to employ the efficient re- 
cursive Green-function method for banded Hamiltonian 
matrices |7j to find the two-particle localization length. 
We obtain the exact Go S, 



(n, n\G$(E)\m, m) 



X, 5—5 . (8) 



ho 



E — Ei — Ej 



by solving the Anderson model in the absence of the in- 
teraction. Here the tfii are the exact single-particle wave 
functions. Clearly, our method of computing £2 is accu- 
rate whenever the enhancement factor £2 /£i is sufficiently 
large. Deviations from the exact £2 arise for small u and 
small £1 where the enhancement is weak. For small u this 
can be easily seen because the second factor in (^) gives 
lim^^o £2 = 0; while the exact limit is £i/2. 

We have studied the two-particle localization length £ 2 
for both fermions and bosons for one-particle localization 
lengths 4.2 < £1 < 105 and interaction strengths < 
u < 1. We find that it is sufficient to use systems with 
N = 500 sites except for the two largest values of £1 where 
we used N = 1000. For each value of £1 we averaged over 
50 realizations of the disorder. In Fig. 1 we have plotted 
£2/^1 as a function of u£%/t at the center of the band 
(E=0) . We have included data for ten values of u for each 
of five values M of £1. The observed scaling behavior 
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is the central result of this paper. While the data in Fig. 1 
are for E — 0, we find that the same scaling behavior 
holds also away from the center of the band |To| . 

This scaling behavior implies that the scale for the in- 
teraction strength u is the energy i/£i which can be in- 
terpreted as the single-particle level spacing = vrt/£i 
of a wire of length £1. The scaling function in Fig. 1 can 
be fit well by a straight line for sufficiently large u/A^ 
while there are deviations from linear behavior for small 
values of the scaling variable. It is natural to suppose 
that these deviations from linear behavior are an artifact 
of disregarding the first factor in Eq. (Q) in computing 
£2- In fact, as shown by the full line in Fig. 1, the linear 
behavior for large u/A^ extrapolates to £2/6- = 1/2 for 
it/Aj — > as expected for the exact two-particle localiza- 
tion length £2 computed from the full expression for G. 
Hence, one may conjecture that the exact scaling relation 
at the center of the band has the form 
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where G ~ 0.17 for bosons and G ~ 0.18 for fermions 
]TT| . We used the fact that at the center of the band there 
is an exact symmetry between attractive and repulsive 
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FIG. 1. Scaling plot £2/^1 = f(u£i/t) for the two-particle 
localization length £2 of (a) bosons and (b) fermions as a 
function of interaction strength u and one-particle localiza- 
tion length £1. Ten values of u are included for each of the 
five values of disorder W = 5 (pluses), W = 4 (squares), 
W — 3 (triangles), W — 2 (crosses), and W = 1.5 (dia- 
monds). The full lines show that the linear behavior for large 
u£i/i extrapolates to £2/^1 = 1/2 for u — *■ 0. The deviation 
from linear behavior for small u^i/t is most likely an artifact 
of our numerical method. 

interactions so that £2 depends only on the absolute value 
of u. Note that while this result is presumably valid for 
arbitrary values of £1, one expects the scaling relation 
to break down for large u where u/t should become a 
relevant parameter due to density-of-states effects. 

Shepelyansky's original prediction, Eq. (|l|), is not con- 
sistent with the scaling relation (|J). E.g., at the cen- 
ter of the band our results show that while £2 depends 
quadratically on £1 as previously predicted, it exhibits 
an unexpected linear (instead of quadratic) dependence 
on the interaction strength u. More generally, the scal- 
ing (||) implies that the enhancement effect sets in for 

weaker interactions u ~ t/£i [compared to u ~ t/^ 2 
according to Eq. ([[])] than originally predicted. This 
result is surprising in view of the following estimate. 
It may be argued [|| that £2 should deviate from £1 
once the two-particle product states \<pi,<f>j) of (unper- 
turbed) energy Eij are strongly mixed by the interac- 
tion U . According to perturbation theory, strong mixing 
occurs when (fa, fa\U\fa, fa)/(Ei : 2 — ^3,4) is of order 
unity. Each \fa,fa) is typically coupled appreciably to 
£1 states. The corresponding matrix element can be es- 
timated as u/^ 2 and the energy denominator as 

Thus, according to this estimate, strong mixing 
1/2 

occurs once it ~ V£i ■ Since this interaction strength 
is large compared to Ag, a comparison with our result 
would suggest that, surprisingly, strong mixing of the 



two-particle product states is not necessary for the en- 
hancement of the two-particle localization length. 

Originally, Shepelyansky [Q approached the problem 
by an approximate mapping to a banded-random-matrix 
model. We have also investigated an alternative random- 
matrix model which is suggested by Eq. (Q) due to 
the band-matrix structure of Go- An extension of 
this banded-random-matrix model will be applied be- 
low to study two interacting particles in a quasi-one- 
dimensional wire. In contrast to ordinary banded random 
matrices we find that the matrix elements g of Gq(E = 0) 
have a Cauchy distribution, P(g) = (r/7r)/(r 2 + g 2 ). To 
obtain this distribution function we argue as follows. For 
defmiteness, consider a diagonal matrix element of Go- 
Due to the localized nature of the wave functions there 
are of the order of £ 2 terms in the sum in (ft). Fur- 
thermore, normalization implies that the wave functions 
are of order l/Vfi within a region of size £1. Hence, 

(n,n\G \n,n) ~ (1/Ci) EEi(V^) ; where x k E.-E, 

is a random variable in the range — At 5; x k ~ At. Thus, 
the matrix elements of the Green function are given by 
averages over random variables whose second moments 
diverge. Neglecting correlations between the Xk, the 
central-limit theorem implies for sufficiently large £1 that 
the diagonal matrix elements have a Cauchy distribution 
of width r ~ 1/t Jl3| . The same argument can be made 
for the off-diagonal matrix elements. Their width T is 
reduced by a factor of exp(— 2\n — m|/£i). As shown in 
Fig. 2, these conclusions arc well supported by numerical 
results. 

We argue that, for sufficiently large u, the resulting 
banded-random-matrix model predicts £2 ~ £1 in agree- 
ment with Eq. (jl]) when neglecting correlations between 
the matrix elements. This result would follow immedi- 
ately from analytical results for banded random matrices 
if the distribution of the matrix elements had a finite vari- 
ance. In this case, the localization length is proportional 
to the square of the bandwidth The same result 

holds true for banded Cauchy matrices for the follow- 
ing reason: Since the eigenstates of the banded Cauchy 
matrix are localized, they effectively sample only a fi- 
nite number of matrix elements drawn from the Cauchy 
distribution. Hence, there exists a corresponding typi- 
cal largest matrix element g max [0|. Beyond g max the 
Cauchy distribution can be cut off, and the resulting ef- 
fective distribution of matrix elements has a finite vari- 
ance. This implies that banded Cauchy matrices belong 
to the same universality class as ordinary banded random 
matrices. We confirmed this conclusion numerically by 
computing the localization length of banded Cauchy ma- 
trices as a function of bandwidth as shown in the upper 
inset of Fig. 2. We have also studied the u dependence of 
£2 predicted by this banded-random-matrix model. How- 
ever, we find that the £2 computed from this model does 
not exhibit the scaling (p|) found for the exact solution. 



3 




e 

FIG. 2. Distributions of diagonal and off-diagonal matrix 
elements g of the projected two-particle Green function Go 
for £1 =46.6 (solid lines). The dashed lines are fits to Cauchy 
distributions. Their width F decreases exponentially away 
from the diagonal on the scale £i/2 as shown in the lower 
inset. The upper inset shows £2 vs £1 as obtained from the 
banded-random-matrix model discussed in the text. These 
data obtained for systems with 10 6 sites confirm our argu- 
ments that £2 ~ £1 in this model. 

Presumably, this is due to correlations between the ma- 
trix elements of Go- For example, correlations in the 
exact Go are implied by Eq. (g) for exceptionally large 
matrix elements. Large matrix elements are due to small 
energy denominators. Each product state with energy 
close to E leads to a correlated £1 x £1 block of large 
matrix elements in Go- Such correlations are neglected 
in the banded-random-matrix model. 

It is an interesting problem to study two-particle lo- 
calization in more than one dimension. In the absence 
of an understanding of the physical origin of the scaling 
parameter u/Ag, it is not clear how to generalize our 
results to these cases. For quasi-one-dimensional wires 
with a finite number of channels M, one easily derives 
a generalized banded-random-matrix-model. Assuming 
that this banded-random-matrix model again correctly 
predicts the dependence of £2 on the bandwidth and com- 
bining the result with scaling suggests that the scaling 
function remains linear. For quasi-one-dimensional wires 
we can order the doubly-occupied sites sequentially along 
the longitudinal direction. When £1 is larger than the 
transverse dimensions of the wire, the bandwidth of Go 
is equal to M£i yielding a corresponding "localization 
length" (M£i) 2 . The actual localization length in the 
longitudinal direction is smaller by a factor M, hence 
£2 ~ M£,i- Finally, assuming the above scaling behavior 
and noting that ~ t/M£i, we obtain a linear scaling 
function, £2/6. ~ M/A^. 



In summary, we have studied the interaction-induced 
delocalization of two particles in a one-dimensional ran- 
dom potential by a novel and efficient numerical ap- 
proach. We have found that the two-particle localiza- 
tion length £2 for coherent propagation of the two par- 
ticles satisfies the scaling relation £2/^1 = /(m/Aj) as 
a function of interaction strength u and one-particle lo- 
calization length £1. This implies that the effect sets 
in for weaker interactions than previously predicted. At 
the center of the band our data suggest that the scal- 
ing function is linear. At present, we do not have a good 
physical understanding of this unexpected scaling behav- 
ior. It will be interesting to see whether the scaling found 
in this paper can be generalized to higher dimensions or 
whether it is a specific feature of one dimension. 

It would also be interesting to study implications of 
coherent propagation due to interactions at finite parti- 
cle density. Our numerical approach can be extended to 
study the propagation of quasiparticlc pairs in the An- 
derson insulator. This will be the subject of a separate 
publication. 

One of us (FvO) is particularly indebted to B. Altshu- 
ler and A. Luther for important discussions which even- 
tually led to a significant revision of a previous version 
of the manuscript. We also acknowledge useful discus- 
sions with G. Hackenbroich, B. Huckestein, A. Muller- 
Groeling, D.L. Shepelyansky, and H.A. Weidenmuller. 



[1] For a review, see B. Kramer and A. MacKinnon, Rep. Prog. 

Phys. 56, 1469 (1993). 
[2] For a review, see D. Belitz and T.R. Kirkpatrick, Rev. 

Mod. Phys. 66, 261 (1994). 
[3] O.N. Dorokhov, Zh. Eksp. Teor. Fiz. 98, 646 (1990) [Sov. 

Phys. JETP 71, 360 (1990)]. 
[4] D.L. Shepelyansky, Phys. Rev. Lett. 73, 2607 (1994); F. 

Borgonovi and D.L. Shepelyansky, preprint (1995). 
[5] Y. Imry, Europhys. Lett. 30, 405 (1995). 
[6] K. Frahm, A. Muller-Groeling, J.-L. Pichard, and D. Wein- 

mann, Europhys. Lett. 31, 169 (1995); D. Weinmann, A. 

Muller-Groeling, J.-L. Pichard, and K. Frahm, preprint 

(1995). 

[7] B. Huckestein, Rev. Mod. Phys. 67, 357 (1995). 

[8] Note that this Green function is properly symmetrized al- 
though this may not be apparent because of the projection 
onto doubly-occupied sites. 

[9] We have not included our data for W = 1 (£1 = 105) 
in Fig. 1. While they are consistent with a linear scaling 
function /, the corresponding error bars are rather large 
because £2 is comparable to the system size. 
[10] F. von Oppen and T. Wettig, unpublished. 
[11] Note that the dependence of £2 on £1 involves both a linear 
and a quadratic term. From a double-logarithmic plot of 
£2 vs £1 over a limited parameter range, this dependence 



4 



was identified in Ref. |6| as a power law with non-integer 
exponent between one and two. 
[12] E.g., see J. -P. Bouchaud and A. Georges, Phys. Rep. 195, 
127 (1990). 

[13] Y.V. Fyodorov and A.D. Mirlin, Phys. Rev. Lett. 67, 2405 
(1991). 



5 



